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Abstract. We construct the biharmonic heat kernel of a suitable closed self- 
adjoint extension of the bi-Laplacian on a manifold with incomplete edge singu- 
larities. We establish mapping properties of the biharmonic heat operator and 
derive short time existence for certain semi-linear equations of fourth order. We 
, apply the analysis to a class of semi-linear partial differential equations in partic- 

ular the Cahn-Hilliard equation and obtain asymptotics of the solutions near the 
£^ , edge singularity. 

c . 

d . l. Introduction 

1—5 . 

The Cahn-Hilliard equation was proposed by Cahn and Hilliard in [Cah61], 
[CaHi58] as a simple model of the phase separation process, where at a fixed tem- 
perature the two components of a binary fluid spontaneously separate and form 
domains that are pure in each component. Let A denote a self-adjoint extension of 
a Laplace operator. Then the Cahn-Hilliard equation may be stated in the following 
form 



Oh 

on 

4 

d t u(t) + A 2 u(t) + A(«(t) - u(t) 3 )) = 0, u(0) = u Q , te (0, T). 

Global existence for solutions to the Cahn-Hilliard equation has been established 
by Elliott and Songmu [EIS086], and Caffarelli and Muler [CaMu95]. In the setup 
of singular manifolds however, there is still a question of asymptotics of solutions 
at the singular strata. This aspect has been studied by the recent work of Roidos 
and Schrohe [RoScl2] in the context of manifolds with isolated conical singularities, 
which has motivated the present discussion here. Using the notion of maximal regu- 
larity, they establish short time existence of solutions to the Cahn-Hilliard equation 
in certain weighted Mellin-Sobolev spaces which then yields regularity and asymp- 
totics of solutions near the conical point. 

In this paper we study the Cahn-Hilliard equation in the geometric setup of spaces 
with incomplete edges, which generalizes the notion of isolated conical singularities. 
Our method is different from [RoScl2] and uses the microlocal construction of the 
heat kernel for the bi-Laplacian. Precise analysis of the asymptotic behaviour of the 
heat kernel allows for a derivation of the mapping properties of the heat operator and 
a subsequent argument on the short time existence for certain semi-linear equations 
of fourth order, in analogy to [JL03] and [BDV11]. Existence and asymptotics 
for solutions to the Cahn-Hilliard equation comes as an application of the general 
results. 

The microlocal analysis of the biharmonic heat kernel on edge spaces also allows 
for derivation of Schauder estimates and ultimately leads to short time existence 
results for fourth order PDE's, including the Bock flow [BaDyll] on singular spaces. 
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This will be the subject of forthcoming analysis. We also point out that our approach 
is not limited to squares of Laplacians on functions, but yields similar results for 
general powers of Hodge Laplacians on differential forms along the same lines. 

Our main result may be summarized briefly as follows. 

Theorem 1.1. Let (M,g) be a Riemannian manifold with an admissible Riemann- 
ian incomplete edge metric g. Then the biharmonic heat kernel H for the Friedrichs 
self-adjoint extension of the Laplace- Beltrami operator A g lifts to a polyhomogeneous 
distribution on a blowup of M + x M 2 and defines a biharmonic heat operator with 
various mapping properties, such that certain semi-linear fourth order equations, 
among then the Cahn-Hilliard equation, admit a short time solution of a certain 
regularity. 

This paper is organized as follows. We first introduce the basic geometry of 
incomplete edge spaces in §2. We then provide a microlocal construction of the heat 
kernel for a certain self-adjoint extension of the bi-Laplacian in §3. We proceed 
with with establishing mapping properties of the biharmonic heat operator in §4. 
We conclude with a short time existence result for certain semi-linear equations of 
fourth order and apply the analysis to the particular example of the Cahn-Hilliard 
equation in §5. 
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2. Spaces with incomplete edge singularities 

We introduce the fundamental geometric aspects of spaces with incomplete edge 
singularities, described in detail in [Maz91], compare also [MaVe12]. Let M be 

o 

a compact stratified space with a single top-dimensional stratum M and a single 
lower dimensional stratum B By the definition of stratified spaces, B is a smooth 
closed manifold. Moreover there is an open neighbourhood U C M of B and a radial 
function x : U D M — > R, such that U (1 M is a smooth fibre bundle over B with 
fibre ^(F) = (0, 1) x F, a finite open cone over a compact smooth manifold F. The 
restriction of x to each fibre is the radial function of that cone. 
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The singular stratum B in M may be resolved and defines a compact manifold M 
with boundary dM, where dM is the total space of a fibration <ft : dM — > B with 
the fibre F. The resolution process is described in detail for instance in [Maz91]. 
The neighborhood U lifts to a collar neighborhood % of the boundary, which is a 
smooth fibration of cylinders [0, 1) x F over B with the radial function x. 

Definition 2.1. A Riemannian manifold (M\dM, g) := (M,g) has an incomplete 
edge singularity at B if over the metric g equals g = go + h, where \h\ go = 0(x) 
as x — > and 

g \ W\dM = dx 2 + x 2 g F + <p*g B , 

with g B being a Riemannian metric on the closed manifold B, and g F a symmetric 
2-tensor on the fibration dM which restricts to a Riemannian metric on each fibre 
F. 

We set m = dimM, b = dimi? and / = dimF. Clearly, m = 1 + 6 + /. We assume 
henceforth / = dimF > 1. Otherwise the incomplete edge manifold reduces to a 
compact manifold with boundary, where our discussion below is no longer applicable. 

Similarly to other discussions in the singular edge setup, see [Alb07], 
[BDVll],[BaVell] and [MaVe12], we additionally require <p : {dM,g F + <p*g B ) 
(B,g B ) to be a Riemannian submersion. If p e dM, then the tangent bundle T p dM 
splits into vertical and horizontal subspaces as dM ®T B dM , where T^dM is the 
tangent space to the fibre of through p and T B dM is the orthogonal complement 
of this subspace. <fi is said to be a Riemannian submersion if the restriction of the 
tensor g F to T B dM vanishes. We say that an incomplete edge metric g is feasible 
if (p is a Riemannian submersion and in addition the Laplacians associated to g F at 
each y 6 B are isospectral. 

The reasons behind the feasibility assumptions are as follows. Let y = (y%, ...,yb) 
be the local coordinates on B lifted to dM and then extended inwards. Let z = 
(zi, Zf) restrict to local coordinates on F along each fibre. Then (x,y, z) are the 
local coordinates on M near the boundary. Consider the Laplace Beltrami operator 
A associated to (M,g) and its normal operator N(x 2 A) yo , defined as the limiting 
operator with respect to the local family of dilatations (x, y, z) — >■ (Xx, X(y — yo), z) 
and acting on functions on the model edge x F x with incomplete edge 
metric g ie = ds 2 + s 2 g F + \du\ 2 . Under the first admissibility assumption, N(x 2 A) yo 
is naturally identified with s 2 times the Laplacian for that model metric. 

The second condition on isospectrality is imposed to ensure polyhomogeneity of 
the associated heat kernels when lifted to the corresponding parabolic blowup space. 
More precisely we only need that the eigenvalues of the Laplacians on fibres are 
constant in a fixed range [0, 1], though we still make the stronger assumption for a 
clear and convenient representation. 

Definition 2.2. Let (M,g) be a Riemannian manifold with a feasible edge metric. 
This metric g = go + h is said to be admissible if in addition to feasibility 

(i) the lowest non-zero eigenvalue Ao > of the Laplace-Beltrami operator 
associated to (F, g F \^( y )) for any y G B, satisfies Ao > dimF. 

(ii) h vanishes to second order at x = 0, i.e. \h\ go = 0(x 2 ) as x — > 0. 
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The reasons behind the admissibility assumptions are of technical rather than geo- 
metric nature and somewhat less straightforward to explain. However we point out 
that condition (i) is satisfied by a rescaling of g F which geometrically corresponds 
to decreasing the cone angle uniformly along the edge. Condition (ii) in particular 
holds for even metrics which depend on x 2 instead of x. The admissibility assump- 
tions yield particular information on the heat kernel expansion, which is then used 
in Proposition 3.5. 

An important ingredient in the analysis of singular edge spaces is the vector space 
V e of edge vector fields smooth in the interior of M and tangent at the boundary DM 
to the fibres of the fibration. This space V e is closed under the ordinary Lie bracket 
of vector fields, hence defines a Lie algebra. Its description in local coordinates is as 
follows. Consider the local coordinates (x, y, z) on M near the boundary. Then the 
edge vector fields V e are locally generated by 

, XOy 1 , . . . , XOy b , Q Z \ ) • • • ) f } ' 

We may now define the Banach space of continuous sections ^°(M), continuous 
on M up to the boundary and fibrewise constant at x — 0. This is precisely the 
space of continuous sections with respect to the topology on M induced by the 
Riemannian metric g. Banach spaces of higher order are defined as follows, compare 
[BDV11]. 

Definition 2.3. Let (M, g) be a Riemannian manifold with an incomplete edge met- 
ric. Let V denote a collection of A and a choice of derivatives in {x~ l V 2 , x _1 V e , V e }, 
which will be specified later. Then for each k G N we define a Banach space 

(M, V) := {u G C 2k (M) n ^°(M) | X o A j u G &?(M),X G V, j = 0, .., k - 1}, 

k 

with the norm ||«||2fc : = || M l|oo + \\X o A- J w|| 00 . 

j=o xev 

3. MlCROLOCAL HEAT KERNEL CONSTRUCTION 

3.1. Self-adjoint extension of the bi-Laplacian. Let A denote the Laplace 
Beltrami operator acting on functions on an incomplete edge space (M, g) with a 
feasible incomplete edge metric g. As a first step we fix a self-adjoint extension of 
the bi-Laplacian A 2 . Consider the space of square-integrable forms L 2 (M,g), with 
respect to g. The maximal and minimal closed extensions of A are defined by the 
domains 

X? max (A) := {u G L 2 (M,g) | Au G L 2 (M,g)}, 

P min (A) := {u G X? max (A) | 3 «j- G C °°(M) such that (3.1) 
Uj u and Auj — > Au both in L 2 (M, g)}. 

where Au G L 2 is a priori understood in the distributional sense. Under the unitary 
rescaling transformation in the singular edge neighborhood 



L 2 (W,dvol(g)) ->L 2 (W,x- f dvol(g)), u^x f/2 u, (3.2) 
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the Laplacian A takes the form 

d 2 1 ( A ff-i^ 2 
1 — ' • 

dx 2 x 2 



A* -foAor 1 -- + — \A Fty + 



Lemma 3.1 ([MaVe12]). Let (M,g) be an incomplete edge space with a feasible 
edge metric. Consider the increasing sequence of eigenvalues (crj)^ =1 of Ap^ y , counted 
with their multiplicities, such that v 2 := o~ + (/ — l) 2 /4 G [0,1). The associated 
indicial roots are given by 7* = ±Uj + 1/2. T/ien any u G £> max (A) expands under 
the rescaling (3.2) as x — > 
p 

$M ~ S (c+[m]t/;+(x, z; y) + C 7[<a]<07(z, y)) + <M, G P min (A), 
i=i 

where the leading order term of each ip^ is the corresponding solution of the indicial 
operator. More precisely, let <pj denote the normalized a j-eigenf unctions of of Ap, y - 
Then 



ij>t{x, z \ y) = xlj (Pj( z i y), i>j {x, z; y) 



j Vx (log x)(j) j (z;y), v$ = 0, 
\ x^ (1 + ajX)4>j(z; y), Vj > 0, 



with cij G R uniquely determined by A. The coefficients Cj[u} are of negative regu- 
larity in y and the asymptotic expansion holds only in a weak sense, i.e. there is an 
expansion of the pairing f B u(x, y, z)x{y) dy for any test function x G C°°(B). 

The Friedrichs self-adjoint extension of A has been identified in [MaVe12] as 

V(A^) = {u G T? max (A) | V j= i,.. iP : cj[u] = 0}. (3.3) 

Note that the sequence of eigenvalues (o~j)? =1 of Ap, y starts with o~j = for j = 
1, .., dim H°(F). The corresponding indicial roots compute to jj = 1/2 ± (/ — l)/2 
and the coefficient are constant in z G F. Consequently, in view of the explicit form 
of the rescaling transformation in (3.2), the Friedrichs domain contains precisely 
those elements in the maximal domain whose leading term in the weak expansion 
as x — > is given by x° with fibrewise constant coefficients. In particular 

X? max (A) n <^°(M) c V(A*). (3.4) 

We fix a self-adjoint extension of the bi-Laplacian as the square of Aj? 

V{A%) = {u G Aj? | Au G A#}. (3.5) 

3.2. Biharmonic heat kernel on a model edge. In the second step we construct 
the heat kernel for A\. We begin with studying the homogeneity properties of the 
heat kernel for the bi-Laplacian in the model case of an exact edge [£ = M. b x 
^(F), dy 2 + g) where i^o(F) = (0, oo) x F, g = ds 2 + s 2 g F ) is an exact unbounded 
cone over a closed Riemannian manifold (F,g F ). The Laplacian A^ on the exact 
edge is then a sum of the Laplacian on (^(F),g) and the Euclidean Laplacian on 
M b . Consider the scaling operation (A > 0) 

^ A : C°°(R + x £ 2 ) -> C°°(R + x £ 2 ), 

(*a«)(*, 0, y, z), (s, y,z)) = u(\H, (As, X(y-y),z), (As, Ay, z)). 
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Under the scaling operation we find 

(d t + A 2 £ )^ x u = \^ x (d t + A 2 £ )u. (3.6) 

Consequently, given the heat kernel H £ for the Friedrichs extension of A| (or at 
that stage any other self-adjoint extension), any multiple of ^\H £ still solves the 
heat equation and also maps into the domain of A|. For the initial condition we 
obtain substituting Y = Xy, S = Xs 

lim / C$ x H £ )(t, s, y, z, s, y, z)u(s, y, z)s f dsdydz 
t^o J £ 

= Ymv\- l - b - f [ H £ (X 4 t, Xs, Xy-Y, z, S, Y,z)u(S/X, Y IX, I)S f dSdYdz 

J £ 

=X- 1 - b - f u(Xs/X, Ay/A, z) = X-^-fuis, y, z). 

By uniqueness of the heat kernel we obtain 

V X H £ = X l - b - f H £ . (3.7) 

In addition to the homogeneity properties of H £ , we also require a full asymptotic 
expansion of the biharmonic heat kernel as {s,s) — > 0. We accomplish this by 
establishing an explicit integral representation of H £ . Under the unitary rescaling 
(3.2) and a spectral decomposition of L 2 (F,g F ) into cx 2 -eigenspaces of A F , we may 
write for the rescaled model edge Laplacian 

a;" -a 




(J 

where is(a) := \/ cr 2 + (f — l) 2 /4 and is defined on C^°(0, oo). The Friedrichs 
extension of A £ is compatible with the decomposition, compare a similar discussion 
in ([Ver08], Proposition 4.9). As a special case of Lemma 3.1, l v has unique self- 
adjoint extension L v in L 2 (IR + ) for v > 1, and in case of v G [0, 1), solutions 
u G ^(^,max) admit a partial asymptotic expansion as s — > 

u(s)=u + c + [u]s^ 1 / 2 + c-[u]\ ' 1 ' ief(U). 

[ Vslog(s), v = 0, 

Then the Friedrichs extension L v of l v is defined, similar to (3.3), by requiring 
c~[u] = 0, and moreover, identifying A R b with its unique self-adjoint extension in 
L 2 (M. b ), we may write 

Af^ = 0A, (<T) +A Ri . (3.8) 

Consequently, it suffices to construct the biharmonic heat kernel for L v + A R & in 
L 2 (IR + x R fe ). Denote by J u the v-th Bessel function of first kind and consider the 
Hankel transform of order v > 



{J^ u u){s):= Vss'J u {ss')u{s')ds', u G C °°(0,oo). (3.9) 
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By ([Co59], Chapter III) and also by ([Les97], Proposition 2.3.4), the Hankel 
transform extends to a self-adjoint isometry on L 2 (IR + ). We denote by 

(^«)(0 := (2tt)- 6 / 2 / u(y)e-*y<dy, ueC™(R b ), 

JR b 

the Fourier transform on R b , which extends to an isometric automorphism of L 2 (R b ). 
Consequently, @ v := Jif u o & defines an isometric automorphism of L 2 (IR + x R b ) 
such that W~ l = Jt^o &~ x . Applying ([Les97], Proposition 2.3.5), we arrive at the 
following 

Proposition 3.2. The isometric automorphism & u diagonalizes L u + A K b. More 
precisely, 

®{L V + A Rb ) = {uG L 2 (R + x R b ) | (S 2 + \E\ 2 ) %u G L 2 (R + x R b )}, 

%(L„ + ArO^T 1 = S 2 + |5| 2 , 

where X, S denote multiplication operators by x G R + and £ G M 6 , respectively. 

Similary, the isometry ^ diagonalizes the squared operator (L u + A R b) 2 , identify- 
ing its action with (S 2 + |S| 2 ) 2 . Consequently we may express the biharmonic heat 
kernel of (L u +A R b) 2 as an integral in terms of Bessel functions. For u G C^°(1R + xR b ) 
we find 

(e-^+^u) (s,y) = (%e-« s2+ \ s V 2 %- l u) (s,y) 

p poo / p poo ^ \ 

= (27r)- 6 / 2 / I e i ^-y^^J v (sp)J v {s P )pe- t{p2+ ^ 2)2 dpdn 

JR b JO \JR b JO J 

u(s, y) d'sdy. 

Denote by <p a the normalized <r 2 -eigenf unction of Ap, where we count the eigenvalues 
cr 2 g Spec(Ai?) with their multiplicities. Then, as a consequence of (3.8), we finally 
obtain for the $-rescaled biharmonic heat kernel on a model edge 

p /*oo 

Hf = {2n)- b ' 2 Q) / ^("^^(^J^^pe^^l 2 ) 2 ^^ 

„ jR b Jo 

■ (j) a (z) <g> <f>a(z). 

The v-th Bessel function of first kind admits an asymptotic expansion for small 
arguments J V (C) ~ Yl'jLo a jC +2 ^ i as C ~~ 0- This yields an asymptotic expansion of 
Hf as (s,s) — > and consequently, rescaling back, we obtain as s — > 

H £ (t, s, y, z, s,y,z)~ ^2a U)j {t,s,y,y,z,z)s 1+2j , (3.10) 

7 

where the summation is over all 7 = — (/ — l)/2 + a/o -2 + (/ — l) 2 /4 with a 2 G 
SpecA^, counted with multiplicity, and each coefficient a v j lies in the corresponding 
cr 2 -eigenspace. We summarize the properties of H?, established above, in a single 
proposition for later reference. 

Proposition 3.3. Consider the model edge (£ = R b x c if(F),dy 2 +g), where (^(F) = 
(0, 00) x F,g = ds 2 + s 2 g F ) is an exact unbounded cone over a closed Riemannian 
manifold [F* , g F ) . Fix the Friedrichs self-adjoint extension of the associated Laplace 
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Beltrami operator As . Then the biharmonic heat kernel He of A| is homogeneous 
of order ( — 1 — b — f) under the scaling operation (X > 0) 

^ A : C°°{R + x S 2 ) -»■ C°°(R + x 8 2 ), 

(*a«)(*, (s, z), (s, y, z)) = u(X 4 t, (Xs, X{y-y),z), (As, Ay, z)). 

Moreover, He admits an asymptotic expansion as (s,s) — > wi/i t/ie mrfez set gwen 
by E + 2N , w/iere 

£ = { 7 > | 7 = -— ^ + V^T^ 2 ' a * G Spec A F }, 

uniformly in other variables and with coefficients taking value in the corresponding 
a 2 -eigenspace. 

3.3. Construction of the biharmonic heat kernel. We can now proceed from 
the analysis of the heat kernel on the model edge to the construction of the heat 
kernel H for the bi-Laplacian on a space (M, g) with an incomplete feasible edge 
metric. The heat kernel construction here follows ad verbatim the discussion in 
[MaVe12] for the edge Laplacian, with the only difference that for the bi-Laplacian 
now rather t 1//4 instead of \ft is treated as a smooth variable. 

The heat kernel is a function on M% = M + x M 2 . Let (x, y, z) and (x, y, z) be the 
coordinates on the two copies of M near the edge. Then the local coordinates near 
the corner in M\ are given by (t, (x, y, z), (x, y, z)). The kernel H(t, (x, y, z), (x, y, z)) 
has a non-uniform behaviour at the submanifolds 

A = {(t = 0, (0,y,z), (0,y,z)) G R + x dM 2 \ y = y}, 

D = {(t = 0, (x, y, z),{x,y,z)) G M + x M 2 \x = x,y = y, z = z}. 

Exactly as in the case of the Hodge Laplacian on edges, see [MaVe12], we introduce 
an appropriate blowup of the heat space M%, such that the corresponding heat kernel 
lifts to a poly homogeneous distribution in the sense of the definition below. This 
procedure has been introduced by Melrose in [Mel93]. For self-containment of the 
paper we repeat the definition of polyhomogeneity as well as the blowup process 
here. 

Definition 3.4. Let W be a manifold with corners, with all boundary faces embed- 
ded, and {(Hi, pi)}f =l an enumeration of its boundaries and the corresponding defin- 
ing functions. For any multi-index b = (&!,..., G we write p b = p^ 1 . . . p b ^ . 
Denote by Vb(2n) the space of smooth vector fields on 2H which lie tangent to all 
boundary faces. A distribution w on QU is said to be conormal, if u G p b L°°(W) for 
some b G and V l . . . V e co G p b L°°(W) for all Vj G H(2U) and for every £ > 0. An 
index set Ei = {(7,p)} C C x N satisfies the following hypotheses: 

(i) Re (7) accumulates only at plus infinity, 

(ii) For each 7 there is P 7 G No, such that (j,p) G Ei if and only if p < P 1 , 

(iii) If (7,p) G E h then (7 + j,p') G E { for all j G N and < p' < p. 

An index family E = (E\, . . . , En) is an A^-tuple of index sets. Finally, we say that 
a conormal distribution u is polyhomogeneous on QU with index family E, we write 
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oj G £^f g (2n), if oj is conormal and if in addition, near each Hi, 

oj ~ ^2 a 7,pPZ( lo gPi) P ' as Pi -> °> 

with coefficients a 7iP conormal on if;, polyhomogeneous with index Ej at any HidHj. 

The homogeneity property (3.7) contains the information how precisely the sub- 
manifolds A, D C Ml need to be blown up such that the heat kernel becomes 
polyhomogeneous. To get the correct blowup of M\ we first bi-parabolically (t 1 / 4 is 
viewed as a coordinate function) blow up the submanifold 

A = {(t, (0,y,z), (0,y,z)) G R + x dM 2 : t = Q,y = y} C M 2 h . 

The resulting heat-space [M%, A] is defined as the union of M%\A with the interior 
spherical normal bundle of A in M\. The blowup [M 2 , A] is endowed with the unique 
minimal differential structure with respect to which smooth functions in the interior 
of M 2 and polar coordinates on M% around A are smooth. As in [MaVe12], this 
blowup introduces four new boundary hypersurfaces; we denote these by ff (the front 
face), rf (the right face), If (the left face) and tf (the temporal face). 

The actual heat-space blowup is obtained by a bi-parabolic blowup of [M 2 , A] 
along the diagonal D, lifted to a submanifold of [M 2 , A]. The resulting blowup 
is defined as before by cutting out the submanifold and replacing it with its spherical 
normal bundle. It is a manifold with boundaries and corners, visualized in Figure 
below. 



£l/4 




Figure 1. Heat-space Blowup 

The projective coordinates on are then given as follows. Near the top corner 
of the front face ff, the projective coordinates are given by 

p = t> , £ = -, £ = -, u = , z, y, z, (3.11) 

P P P 

where in these coordinates p, £, £ are the defining functions of the boundary faces ff, 
rf and If respectively. For the bottom corner of the front face near the right hand 
side projective coordinates are given by 

* x y-y — ,q n\ 

r = — , s = -, u = — , z, x, y, z, (3.12) 

X 4 X X 

where in these coordinates 7", Sj X 3X6 the defining functions of tf, rf and ff respectively. 
For the bottom corner of the front face near the left hand side projective coordinates 
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are obtained by interchanging the roles of x and x. Projective coordinates on ^ 2 
near temporal diagonal are given by 

t 1 / 4 (x-x) rj_y-y 7 _ x{z-i) 

In these coordinates tf is the face in the limit 1(5, U, Z) \ — > oo, ff and td are defined 
by x, i], respectively The blowdown map (3 : — > M 2 is in local coordinates 
simply the coordinate change back to (t, (x, y, z), (x, y, T)). 

In case the edge manifold is an exact edge (S = M. b x c ^(F),dy 2 + g) where 
(tf(F) = (0, oo) x F f ,g = ds 2 + s 2 g F ), Proposition 3.3 implies that H £ lifts 
to a polyhomogeneous conormal distribution on the heat space blowup, of order 
(— m),m = 1 + b + /, at the front and the temporal diagonal faces, vanishing to 
infinite order at tf, and with the index set at rf and If given by E + 2Nq, where 



E = { 7 > | 7 = -^-y^ + y (/ 4 1)2 + a2 ' ff2 G Spec A^}. 

In the general case of a feasible edge space (M,g), Hs is only an initial parametrix 
and solves the heat equation only to first order. Repeating almost ad verbatim the 
heat kernel construction in case of the edge Laplacian in [MaVe12], we arrive at 
the following 

Proposition 3.5. Let (M m ,g) be an incomplete edge space with a feasible edge 
metric g. Then the lift f3*H is polyhomogeneous on ^ 2 of order (— dimM) at ff 
and td, vanishing ot infinite order at tf, and with the index set at rf and If given by 
E + Nq where 



E = {7 > | 7 = "^-y^ + V ^-y^ + cr 2 , \ E Spec A^}. 
More precisely, if s denotes the boundary defining function of rf, we obtain 

(oo oo \ 

s^a^{P*H) + s~< +2+3 a'^{P*H) as s 0, 
3=0 j=0 J 

where the coefficients a 1 j(H) are of order (— m) at the front face and lie in their cor- 
responding Ap y eigenspaces. The higher coefficients a'^ j((3*H) are of order (— m+1) 
"1 .11 

Proof. Recall the heat kernel construction in [MaVe12], which we basically follow 
here. Denote by A the Laplace Beltrami operator on (M, g). We write £, := dt + A 2 
for the heat operator. The restriction of the lift (3*(tC) to ff is called the normal 
operator Ns(tC) y() at the front face (at the fibre over y G B) and is given in 
projective coordinates (3.12) explicitly as follows 

N„(tc) m =tU+ (-% - fs-% + s- 2 A F , m + &fy\ 
=:,(a r+ (A2:> + Ar) 2 ). 

Nff(tC) does not involve derivatives with respect to (yo,x,y,z) and hence acts tan- 
gentially to the fibres of the front face. Consequently in our choice of an initial 
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parametrix Hq we note that the equation 

N s {tC) o N g (H ) = 

is the heat equation for the bi-Laplace operator on the model edge ^(F) x R 6 . 
Consequently, the initial parametrix Hq is defined by choosing Nq(Hq) to equal the 
fundamental solution for the heat operator N s (t£), and extending N S (H ) trivially 
to a neighbourhood of the front face. More precisely, consider the biharmonic heat 
kernel H gy on the model edge {^{F) x R 6 , ds 2 + s 2 gy {) + du 2 with the parameter 
2/o G B. Then in projective coordinates (r, s, yo, z, x, u, 2) near the right corner of ff, 
where x is the defining function of the front face, we set 

H (r,s,u,y ,z,z) := x~ m <f)(x)H gi y (T, s,u, z,s = l,u = 0,z), (3.14) 

where is a smooth cutoff function, = 1 in an open neighborhood of x — 0, and 
with compact support in [0,1). By Proposition 3.3, our initial parametrix Hq is 
poly homogeneous on M 2 and solves the heat equation to first order at the front face 
ff of Moreover Proposition 3.3 asserts 

oo 

H * ~ H J2 ^ +2j ^j( H o), s 0, (3.15) 

7G-B 3=0 

with each coefficient a 7 j(i?o) lying i n the corresponding A F>yo eigenspace. The error 
of the initial parametrix H is given by 

F(t£)H = (/3*K) - r(Aj(f + Af) 2 ) H =: P . 

The leading order term in the expansion of (3*(tC) at td does not depend on the edge 
geometry and corresponds to the bi-Laplacian on a closed manifold. Consequently, 
classical arguments allow to refine the initial parametrix such that the error term 
Pq is vanishing to infinite order at td, compare the corresponding discussion in 
([MaVe12], Section 3.2). We need to understand the explicit structure of the 
asymptotic expansion of Pq at ff and rf. By Definition 2.2 (ii) we find 

/T A 9 = x- 2 (AJJ + Af) + 5-% + L 2 , (3.16) 

where L\ is comprised of the derivatives d u d z and L 2 consists of edge derivatives V 2 . 
Consequently, we obtain after taking squares 

F{tAl) - r(A5f + Af) 2 = rx ((Ajjf + Af )L, + L X (A^ + Af)) 

+ rx 2 ((AJ + Af )L 2 + L 2 (Aj(f + Af ) + Ll) + rx^L\. 

We now apply each of the summands above to the asymptotic expansion (3.15) of 
Hq. Note that A s ,y annihilates each s 7 a 7i o(-ffo), 7 G E, and lowers the s-order of 
s 7+2j a^j(H ) by 2, if j > 1. Consequently, we obtain as s — > 

oo 

P = f3*{t£)H ~ x- m+1 sJ+J ' 2c ^r 

7 G£ 3=0 

The next step in the construction of the heat kernel involves adding a kernel H 
to H , such that the new error term is vanishing to infinite order at rf. In order 
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to eliminate the term s k a in the asymptotic expansion of Pq at rf, we only need to 
solve 

{-d 2 s - fs'% + s- 2 A F<yo ) 2 u = s k {r~ l a). (3.17) 

This is because all other terms in the expansion of tC at rf lower the exponent 
in s by at most one, while the indicial part lowers the exponent by two. The 
variables (r, u, x, y , y, z) enter the equation only as parameters. The equation is 
solved by Mellin transform as well as spectral decomposition over F. The solution 
is polyhomogeneous in all variables, including parameters and is of leading order 
(k + 4). Consequently, the correcting kernel H' Q must be of leading order 2 at rf and 
of leading order (— m + 1) at ff, since Pq is of order (—2) at rf and (— m + 1) at ff and 
the defining function x of the front face enters (3.17) only as a parameter. Hence 

too oo \ 

£ s^a^H,) + £ s^a'^(Hi) as s -> 0, 
3=0 j=0 J 

where the coefficients a 7 j(Hi) each lie in the corresponding A F , yo eigenspace. In the 
following correction steps the exact heat kernel is obtained from Hi by an iterative 
correction procedure, adding terms of the form Hi o (Pi) k , where Pi := tCHi is 
vanishing to infinite order at rf and td. This leads to an expansion 



DC 



0*H ~ [J2 sJ+2j(1 iAP* H ) + Yl s 7+2+ H,i(£*#) as s -> 0, (3.18) 

~(£E \j=0 j=0 / 

where the coefficients a 7J (if) still lie in their corresponding A. F>yo eigenspaces, and 
are of order (— m) at the front face. The higher coefficients a' 1: j((3*H) are of order 
(-m + 1) at ff. 

Note that in the construction above, we have only used feasibility of the edge 
metric and Definition 2.2 (ii). The assumption of Definition 2.2 (i) for admissible 
edge metrics is not required there, but plays an essential role in the argument that H 
indeed takes values in ^(A^). First note that H indeed takes values in f^(A ^), since 
the expansion (3.18) satisfies the characterization of maximal solutions in Lemma 
3.1 and also the condition in (3.3) under the rescaling $. 

The conclusion that AH also takes values in i^(Aj?) is more intricate. Recall 

(3.16). It is easily checked from (3.18) that x~ 2 (Af,^ } + Af ) H indeed takes 

values in @(Ajr) without any further assumptions. Application of (x^Li + L 2 ) to 
H preserves the expansion (3.18), however the coefficients in the expansion need not 
lie in the correct A F ^ yo eigenspaces, and hence we cannot deduce that (x~ 1 Li + L 2 )H 
maps into i^(Aj?) in general. 

By condition (i) of Definition 2.2, any 7 7^ is automatically 7 > 1, and hence it 
then suffices to check whether the s° coefficient in the expansion of (x _1 Li + L2)H 
lies in the zero-eigenspace of A F , in other words is harmonic on fibres and hence 
constant in z. The leading term s ao } o{f3*H) in the expansion of /3*H is annihilated 
by {sd s ), d z and is increased by (3*xd y = sd u + xsd y . Consequently, (x~ 1 Li + L 2 )H 
admits no s° term and hence trivially maps into f^(Ajr). 
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The kernels H and AH thus both map into ^(Ajr) and hence by definition, 
H indeed maps into f^(Ajf) and thus is the biharmonic heat kernel associated to 
A%. ' ' □ 

4. Mapping properties of the biharmonic heat operator 

In this section we prove boundedness and strong continuity of the biharmonic 
heat operator with respect to Banach spaces introduced in Definition 2.3. 

Theorem 4.1. Let (M m ,g) be an incomplete edge space with a feasible edge metric 
g. Fix the Friedrichs extension Ajr of the corresponding Laplace Beltrami operator. 
Put V = {A} and V = {A, x~ l V c 2 , x~ x V e , V e }, where V e C V e consists locally of all 
edge vector fields where xd y is weighted with functions that are fibrewise constant. 
Then the associated biharmonic heat operator e~ tA & is a bounded map between the 
(weighted) Banach spaces 

e -tA>„ . ^fc(M,P ) -> r^f +1) (m,p). 

Proof. First we prove the statement for k — 0. The explicit structure of the heat 
kernel expansion in Proposition 3.5 implies that for any X £ D applied to the 
biharmonic heat kernel H, the lift (3*(XH) admits the following behaviour near the 
front face of the heat space 

0*(XH) = O ((p rf pif) (p ff p td )- m - 2 PtT) , (4-1) 
where p* denotes a defining function of a boundary face *, * £ {rf, If, ff, td, tf}. 

Consider the lift of the volume form in the various projective coordinates near ff. 
We explify the transformation rules for the volume form near the lower left, lower 
right and the top corner of the front face 

, r, t x y-y 
near left corner: r = — = p t f, s = — = pif, u = , x = pg, y, z, z, 

x 4 x x 

f3*(x* dx &vo\q M (x)) = h ■ x m s* ds du dz, 

~ t „ x „ y — y _ 

near right corner: r = — = p tf , s = - = p rf , u = - , z, x = p ff , y, z, 

JU JU JU 

(3*(x,f dx 6.wo\qui{x)) = h ■ x m r _1 dr dudz, 

near top corner: p = t 1 / 4 = p$, £ = — = p r f, ^ = — = pif, u = - — — , y, z, z, 

P P P 

/3*(x^ dx (1vo\bm(x)) = h ■ p m d^ du dz. 

(4.2) 

The projective coordinates and the transformation rule for the volume form where 
the front and the temporal diagonal faces meet, is as follows 

t 1 / 4 x — x y-y x(z-z) 

V= — = Ptd, S = - rjr , U = ^ T?r , Z = tl/4 , x = p s , y, z, (4 3) 

P*{x f dxdvol dM {x)) = h ■ x m r] m dSdUdZ. 

When we combine the asymptotics of (3*(XH) in (4.1) with the behaviour of the 
volume form in the various projective coordinates (4.2) and (4.3), we find that 
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(3*(XHxf dx dvolgM^)) has a singular behaviour of (pgptd) 2 < ct Conse- 
quently, we may estimate for X E T> and any continuous function, in particular 
for any u E ^(M) 

for some constant C > independent of u. Furthermore, by Proposition 3.5, 
f3*XH ~ aoPrf, as p r f — > for X E D, where ao is fibrewise constant, i.e. independent 
of z E F. Here the fact that for X e ^ _1 V e its d y component is weighted with a fi- 
brewise constant function, is essential. Hence, indeed Xe~ tA &u E t~ 1//4 ^(M). This 
proves the statement for k = 0. The general statement follows from the fact that due 
to (3.4), %f(M,V) C and for any u E A k e~ tA ^u = e - tA ^A k u, by 

uniqueness of solutions to the biharmonic heat equation with fixed initial condition. 

Finally note that while t 1 ^ 4 Xe~ tA ^u is continuous even for X E x^ 1 V e , it need 
not remain fibrewise constant at x = since in general X may include vector fields 
weighted with ^-dependent functions. Hence x~ x V e is replaced by x x V e in T>, where 
Vg C V c consists locally of linear combinations of {xd x ,xd y ,d z }, weighted with 
functions that are fibrewise constant. □ 

Remark 4.2. A particular property of u E ^(M, T>) with {d x , d y , x~ l d z } C T> is 
worth noticing. In the singular neighborhood of the edge, the distance defined by 
the Riemannian edge metric g is equivalent to 

d{(x,y,z),(x,y,z)) = (\x-x\ 2 + \y-y\ 2 + (x + xf\z - z\ 2 ) 1/2 . 
In local coordinates near the edge we find 

u(x,y,z) -u(x,y,z) 
= u(x, y, z) - u(x, y, z) + u(x, y, z) - u(x, y, z) + u(x, y, z) - u(x, y, z) 
= d x u(£, y, z) (x-x) + d y u(x, f y,z)(y-y) + x~ l d z ii{x, y,()x(z-z). 
Consequently we obtain 

\u(x,y, z) — u(x,y,z)\ < \\u\\ 2 (\x — x\ + \y — y\ + (x + x)\z — z\) 

< V2 \\u\\ 2 d((x,y,z), (x,y,z)). 
In other words, u E ^(M,V) is automatically Lipschitz with respect to d. 

Our argument on existence of solutions to certain semi-linear fourth order equa- 
tions in the section below crucially depends on the strong continuity property of 
the biharmonic heat operator with respect to the Banach space ^g k (M, T>). This is 
the content of the next theorem. Note that for strong continuity we will choose a 
different space T>' of allowable operators, smaller than in Theorem 4.1. Beforehand 
we note the following well-known function analytic result. 

Lemma 4.3. Let D be a self-adjoint non-negative unbounded operator in a Hilbert 
space H . Then the following is true. 

(i) The continuously differentiable solution to (dt + D 2 )u = 0, with u(t) E 
3i(D 2 ) for t > and lim^o u(t) = u E H is unique, for any given u E H. 

(ii) For any uq E &(D), we have De~ tD uq = e~ tD Duq. 
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Proof, (i) For s G (0, t] we compute 

d s e~ {t - s)D ' 2 u(s) = -d t e- {t - s)D2 u{s) + e~ {t - s)D2 d s u(s) 

= e-^ D2 D 2 u{s) - e-^ D2 D 2 u(s) = 0. 

Consequently, e~( t ~ s ^ L)2 u(s) is constant for s G (0, t]. Since e~ tr>2 converges to Id in 
the Hilbert space norm as t — > and, by assumption, u(t) is continuous at t — 0, we 
find that e~^~ s ^ D u(s) is constant for s G [0, t). Considering the limit of e-^- s ) D u(s) 
as s — > t and as s — > proves for any uq G H 

u{t) = e~ w2 u . 

(ii) Consider A G Res(Z) 2 ) in the resolvent set of D 2 . Then for any uq G S'(D), the 
resolvent (D 2 — A) _1 -u G £)(D 2 ) and we compute 

(-D 2 - A) J D( J D 2 - A)- x wo =.D( J D 2 - X)(D 2 - A)" 1 ^ = .Dmo, 

^D(D 2 - A) -1 m = - A) _1 Du - 

The statement now follows by closedness of D and definition of the heat operator 
as the strong limit e~ tD := lim^oo (J + tD 2 /n) n . □ 

Theorem 4.4. Let (M m ,g) be an incomplete edge space with a feasible edge metric 
g. Fix the Friedrichs extension Aj? of the corresponding Laplace- Beltrami operator. 
Put V = {A,V 2 ,V C }. Then the associated biharmonic heat operator e~ tA & is a 
strongly continuous bounded map between Banach spaces 

Proof. By (3.4), we have tf? k {M,V') C @(A%) and hence for any u G <tf 2k (M,V') 
we infer by the previous Lemma 4.3, A k e~ tA -^u = e~ &A k u. This reduces the 
statement to k — 1 and = 0. Proof of both cases requires stochastic completeness 
of the biharmonic heat kernel, which we explain below. Solutions to the initial value 
problem 

d t u + A 2 u = 0, u(0) = m , u(t) G @{A%), t > 0, 

are unique and in fact given by u(t) = e~ tA &u G 3>(A 2 ^). We have observed 
in subsection 3.1 that, reversing eventual rescaling, the Friedrichs domain contains 
precisely those elements in the maximal domain whose leading term in the weak 
expansion as x — > is given by a; , with a fibrewise constant coefficient, cf. (3.4). 
Consequently, the constant function 1 G St(A^). Moreover, Al = G @{A&) and 
consequently 1 G ^(A 2 ^). The constant function 1 solves the heat equation and 
hence we deduce by uniqueness of solutions the stochastic completeness 

e -t±% i = f H{t, p, p) dvol g (p) = 1, for all p G M, t > 0. (4.4) 
Jm 

This reduces the case to k = 0, 1. We can now prove the statement for k = 0, 
basically repeating the arguments in ([BDV11]) where the classical proof of strong 
continuity of the heat operator on closed (non-singular) manifolds is adapted to the 
present setup. Using stochastic completeness we find 

(e~ tA2 ?u)(p,t) — u(p) = / H (t,p,p) (u(p) — u(p)) dvol g (p). 

Jm 
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Consider the distance function d(p,p) induced by the incomplete edge metric g. In 
the singular neighborhood of the edge, the distance is equivalent to 

d((x,y,z), (x,y,z)) = (\x-x\ 2 + \y-y\ 2 + (x + x) 2 \z - z\ 2 ) 1/2 . 

Note that u G ^(M) is continuous with respect to the topology induced by the 
Riemannian metric g and hence by the distance function d. Hence for any e > 
there exists some 5(e) > 0, such that for d(p,p) < 5(e) one has \u(p) — u(p)\ < e. 
For any given e > we separate the integration region into 

M+ := {p\d(p,p)>5(e)}, 

M~ := {p\d(p,p) <5(e)}. (4.5) 
Employing continuity of u we find 



i -tA 2 



u\ 



H (t,p,p) (u(p) — u(p)) dvol g 



M 



< 



\H (t,p,p) | • \u(p) — u(p)\ dvol g (p) + / \H (t, p,p) | • \u(p) — u(p)\ dvol 



< 2 



£l/4 

5(e) 



\ u \\o 



\H(t,p,p) 



i^7^ 



d(p,p) 

f l/4 



dvoL 



+ e 



\H (t,p,p) | dvolg 



It may be checked in the various projective coordinates around the front face in the 
heat space that (3*(\H\ dvol g ) and (3*(d(p,p)t~ 1 ^)p ti is bounded. Since f3*\H\ 
is vanishing to infinite order at tf, we find that both integrals above are bounded 
uniformly in (t,p,e). Therefore we obtain 



-tA 



u — u\\ < C 



f l/4 

5(e) 



\u\\ + e 



|e tA u — u\\ 



< 2eC for t 1 / 4 < e5(e) 



Mlo- 



Thus, for any given e > we can estimate 
This proves strong continuity of the biharmonic heat operator on ^(M). It remains 
to prove the case k — 1. Strong continuity of the biharmonic heat operator with 
respect to ^ 2 (M) means \\X(e~ tA ^u — u)\\ — > as t — > 0, for u G ^ 2 e (M) and 
X G T>'. If X = A, this follows from the case k = 0, since Ae~ tA ^u = e~ tA ^Au for 
u G ^e(M) C ^(A^). For X G {V e 2 , V e } the leading order of (3*H at the front face 
is preserved under X, so that away from td, the estimates reduce to the case k = 0. 



Near td, a priori X leads to p^ 2 singular behaviour. However, integration by 
parts, exactly as worked out in detail in [BDV11] allows to pass X to u, so that 
the estimates again reduce to the case k = 0. We write down the argument for 
completeness. The edge vector fields obey the following transformation rules in 
projective coordinates (3.13) near the temporal diagonal 

P*(xd x ) = -rjd v + -d s + Zd z + xd x , f3*(xd y ) = -du + xd y , /3*(d z ) = ~d z + d z . 

rj i] 1] 

By Proposition 3.5 

0*H( V , S, U, Z, x, y, z) = x" m r / - m G(r / , S, U, Z, x, y, z), 

(3*(x f dxdvol dM (x)) = h(x7]) m (l - 7]S) f dSdUdZ, 
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where G is bounded in its entries, and in fact infinitely vanishing as \ (S, U, Z)\ — > oo, 
and h — h(r), x(l — r)S),y — xrjU, z — rjZ, x, y, z) is a bounded distribution on 
We consider \\xd x (e~ tA u — u)\\ . Using stochastic completeness of the heat kernel, 



we find 



F : xS . A e^u - u) = /( lS , B )« fifl i^dvoWi) 
(xd x )[Hu(x, y, z)x? dx (1voIqm(x)] ='■ Fi — F 2 . 



Next we transform to projective coordinates and integrate by parts in S, where the 
boundary terms lie away from the diagonal and hence are vanishing to infinite order 
for t — > by the asymptotic behaviour of the heat kernel. Omitting these irrelevant 
terms, we obtain 



Fi = J (-<nd ri + Us + Zd z + xd^j [(xr])- m G(r],S,U,Z,x,y,z)] 
x u(x(l - r]S),y - xrjU, z - rjZ) h(xr]) m (l - r]S) f dS dU dZ 
= J [(-vd v + Zd z + xd x )(xr))- m G] ■ uh(xr)) m (l - r)S) f dS dU dZ 



-d s ) a 

V 



J {xr])- m G ■ u 



h(l-r]S) f dSdUdZ 

-d s ) h(xrj) m (l-rjSy 
V J 



dSdUdZ. 



We perform similar computations for F^. 



(xd x H)u(x, y, z) + H(xd x u) x-^dx dvol^M^ 



-rjd v + -d s + Zd z + xd x 

7] 



(xr])- m Gj u ■ h(xr]) m {l - r]S) f dSdUdZ 
+ J G(r],S,U,Z,x,y,z)(xd x u(x,y,z))h(l-r]S) f dSdUdZ 
= J [(-vdr, + Zdz + xd x )(xri)~ m G]-uh(xri) m (l-riSydSdUdZ 
-J(xr]y m G-u Qd s ^ h(xr]) m (l -r]S) f 



dSdUdZ 



+ J G(r), S,U, Z,x,y, z)(xd x u(x,y, z))h{\ — i]S)^dS dU dZ. 
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Thus F = F\ — F2 becomes 



F = J [(-7 1 d r] + Zd z + xd x )(xr ] y m G(r ] ,S,U,Z,x,y,z)jh(x7 1 ) m (l-7 1 S) 
x (u(x(l — r)S),y — xr]U, z — rjZ) — u(x, y, z))dS dU dZ 



- / G(r],S,U,Z,x,y,z) 



x (u(x(l — T]S),y — xt]U, z 



G 



-d s u(x(l -r]S),y 



-d s )h-(l-r ] Sy 
V J 

- r/Z) — u(x, y, z))dS dU dZ 

— xrjU, z — r]Z) + xd x u(x, y, z) 



h{l- v S) f dSdUdZ. 



Now, each of the three integrals is estimated as above for k = by separating the 
integration region into M e + and M~ for any given e > 0. Note that in the final 
integral we use the fact that u G ^(M, V) so that r/~ 1 dsu and xd x u are bounded. 
Higher order and other edge derivatives may be estimated in a similar way. This 
proves strong continuity in general and as a trivial consequence boundedness of the 
biharmonic heat operator. □ 

Remark 4.5. We point out that Theorem 4.1 holds also when the basic space 
^g(M) is replaced by the Banach space of sections continuous up to x = 0, 
without the requirement of being fibrewise constant at DM. Also, we may set 
V = {A, x _1 Vg , x~ l V e , V e }. The use of the refined space %®(M) and the restriction 
of X" 1 V c to x- x V' e in V becomes however crucial in Theorem 4.4. 

5. Short time existence of semi-linear equations of fourth order 

In this section we explain how the mapping properties of the biharmonic heat 
operator and its strong continuity yields short-time existence of solutions to certain 
semilinear equations of fourth order. The underlying idea is based on the following 
theorem. 

Theorem 5.1. [Tay96, Proposition 15.1.1] Let P be some, possibly unbounded, 
linear operator in a Hilbert space H, bounded from below. Suppose that V,W C H 
are Banach spaces, such that P : V —tW is bounded and moreover 

(i) e~ tp : V — > V is a strongly continuous semigroup, fort > 0. 

(ii) Q : V — > W is locally Lipschitz, 



111 



-tp 



: W — >■ t bounded for some 7 < 1. 



Then for any u e W , the initial value problem 

8 t u -Pu = Q(u), n(0) = u e W 

has a unique solution u 6 C([0, T], ^} k+1 \M, T>)), for some T > 0, where T may 
be estimated from below in terms of \\uo\\w The solution u is the fixed point of the 
operator F : V — > V with 



F(u) 



e~ tp u + 



-(t-s)P 



Q(u)ds. 



Corollary 5.2. Let (M,g) be an incomplete edge space with a feasible edge metric 
g. Put V = {A,V e 2 ,V c j and V = {A, x~ l V,?, x~ x V e , V e }, where V' e C V e consists 
locally of linear combinations of {xd x , xd y , d z }, where xd y is weighted with functions 
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that are fibrewise constant. Assume Q : fe+1 ^ (M, V) — > ^ k (M,V) is locally 
Lipschitz. Then the semilinear equation 

d t u + A 2 u = Q(u), u{0) = u e (M, V) 

has a unique solution u G C([0,T],^ c ( ' k+1 \M,V)), for some T > 0, where T may 
be estimated from below in terms of ||wo||2fc- 

Proof. Consider first a slightly smaller set of operators V = {A,Vg,V e } and set 
W = <g? k (M, V'),V = ^ (fc+1) (M, V). In view of Theorem 4.1 and Theorem 4.4, the 
heat operator associated to A|? satisfies the conditions of Theorem 5.1 with 7 = 1/4. 
Consequently, by Theorem 5.1 the unique solution u exists and lies in ^^ k+1 \M, V). 

This solution is the fixed point of the map F : ^ (fe+1) (M, V) ->• ^ (fc+1) (M, V) 
with ^ 

F(u) = e~ tA ^u + [ e- {t ~ s)A ^Q(u)ds. 
Jo 

However, by Theorem 4.1, F actually maps ^ (fc+1) (M, £>') C (M,V ) to 

^ e {h+l) (M,V). Consequently, u G ^} k+1) (M,V), with a slightly better regularity, 
as claimed. □ 

We now apply this general existence result to the example of the Cahn-Hilliard 
equation on an incomplete edge manifold. We define 

Q : ^f +2 (M,P) -> ^(M,V), Q(u) := A(u - u 3 ). 
The mapping Q is indeed locally Lipschitz, since for any u, v G ^g k+2 (M, X>) 
||Q(u-u)|| 2 fc < ||A(u-i;)||2fc+ ||A(M-w) 3 || 2fc 

< ||« - ^Ikfc+i) (! + Ik - ^IIW)) • 

We hence arrive at our final result. 

Corollary 5.3. Let (M,g) be an incomplete edge space with a feasible edge metric 
g. Put V = {A,V 2 ,V e } and V = {A, ar 1 V 2 , x~ x V e , V e }, where V e C V e conszsfo 
locally of linear combinations of {xd x , xd y , d z }, where xd y is weighted with functions 
that are fibrewise constant. Then the Cahn-Hilliard equation 

d t u + A 2 u + A(u - u 3 ) = 0, u(0) = u G (M, V) 

has a unique solution u G C([0, T], ^ 2k+2 (M, £>)), for some T > 0. 

It should be noted that in correspondence with [RoScl2] our approach leads to 
an explicit identification of the asymptotics of the Cahn-Hilliard solution at x — 0. 
Indeed, u G ^ 2k+2 (M,X>) C ^(A^ -1 ), which yields a partial asymptotics of u to 
higher and higher order, depending on k G N, by an iterative application of Lemma 
3.1 for k steps. 
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